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Abstract 



We show that QCD undergoes dimensional reduction at high temperatures 
also in the quark sector. In the kinematic region relevant to screening physics, 
where the lowest Matsubara modes are close to their "mass-shells" , all static 
Green's functions involving both quarks and gluons, are reproducible in the 
high-T limit by a renormalizable three dimensional Lagrangian up to order 
g 2 (T) ~ 1/lnT. This three dimensional theory only contains explicitly the 
lightest bosonic and fermionic Matsubara modes, while the heavier modes 
correct the tree- level couplings and generate extra local vertices. We also find 
that the quark degrees of freedom that have been retained in the reduced 
theory are nonrelativistic in the high-T limit. We then improve our result to 
order g 4 (T) through an explicit nonrelativistic expansion, in the spirit of the 
heavy quark effective theory. This effective theory is relevant for studying 
QCD screening phenomena with observables made from quarks, e.g. mesonic 
and baryonic currents, already at temperatures not much higher than the 
chiral transition temperature T c . 
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I. INTRODUCTION 



The behavior of field theories, quantum chromodynamics (QCD) in particular, at fi- 
nite temperature (T) is of great phenomeno logical and theoretical interest. In general, we 
can roughly classify finite temperature physics into two categories: real-time dynamics and 
screening phenomena. On one hand, real-time dynamics describes the time- dependent re- 
sponse of a system (correlators as functions of real frequency) to time-dependent external 
probes. On the other hand, screening phenomena refer to the static spatial-dependent re- 
sponse (correlators as functions of the spatial momentum) to time-independent external 
probes. 

Without doubts, real-time dynamics is highly interesting and important, and it could 
apply, in principle, to a wide range of phenomenological applications. Unfortunately, its 
complexity and the present lack of systematic nonperturbative approaches make real-time 
dynamics at finite temperature a territory of field theory that is still, to a very large extent, 
barely cultivated [|TJ. In addition, the difficulties of realizing an equilibrated experiment 
with a given temperature, at least in the specific case of QCD, has practically prevented 
us from accessing real-time data that are not tempered by some ad hoc phenomenological 
assumptions. 

In contrast, the Euclidean nature of the static correlation functions makes the physics 
associated with screening phenomena relatively simpler. In fact, static correlation functions 
are determined from equilibrium ensembles and involve no tricky analytic continuation 
and, therefore, screening physics is well-suited to the lattice formulation of field theories 
at finite temperature. As a consequence, lattice QCD at finite temperature provides us 
with a large body of measurements not only of bulk quantities, such as the specific heat 
and pressure, but also of more detailed observables ranging from screening masses || to 
screening wave functions Q. 

In the past few years several physical pictures or scenarios have been proposed for prop- 
erly understanding and interpreting the available lattice data, both the data involving ob- 
servables made from pure gluonic fields || and the data involving observables made from 
explicit quark fields ||. One of the most important concepts used in these works is the 
so-called dimensional reduction (DR) at high temperatures: this concept can be roughly 
summarized by saying that these QCD screening observables can be effectively described by 
a three-dimensional theory when temperature is high enough. 

The DR picture is based on the observation that two different scales appear, in 
general, in field theories at high T: one scale is order T and the other is order one relative to T 
(this second scale is strictly order one, i.e. independent of T, only at the tree- level in general). 
At high T these two scales become vastly separated and certain degrees of freedom effectively 
decouple; this phenomenon is analogous to the decoupling of heavy particles 0. So far 
the existing literature has mostly concentrated on situations where the observables 

involved are purely bosonic: the zero Matsubara modes (scale of order one) are the explicit 
light degrees of freedom, while the non-zero modes (scale of order T) play the role of the 
heavy degrees of freedom. 

However, we are often interested in observables that couple directly to quark degrees 
of freedom, observables that would vanish without the explicit presence of quarks in the 
theory. Our interest in this kind of observables is by no means academic, but rather it is 
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dictated by the fact that many important observables fall into this class, such as mesonic 
and baryonic correlators, and in fact such observables have been extensively studied at 
T = However, less attention has been paid up to now to such correlators between 
observables made from quark fields in the context of DR. The reason is probably that the 
scale separation in these situations is less clear, since all fermionic Matsubara modes have 
energies of order T, due to their antiperiodic boundary condition. Therefore, only the 
specific underlying dynamics can make DR possible for observables that are made directly 
with quarks. In other words, it is the theory itself that must generate the scale separation, 
a necessary condition for the decoupling of some degrees of freedom. 

In a recent letter |L4] we were able to show in a specific asymptotically free theory, the 
Gross- Neveu model, that a new scale smaller than T, i.e. T/lnT, naturally emerges at 
high T. It is this dynamically generated scale that provides the scale hierarchy that in turn 
makes DR possible for fermionic observables in that model. The purpose of this paper is to 
examine what happens to QCD at high temperature, to verify that QCD undergoes DR to 
one-loop order also in the quark sector, and to give the corresponding explicit form of the 
reduced theory. 

As we will show, the DR in the quark sector is not based on the fact that the lightest 
Matsubara mode (vrT) is much lighter than the next lightest mode (3nT). Instead, the small 
expansion parameter is the so-called off-shellness or residual momentum after subtracting 
7rT. This expansion parameter roughly measures how far is a given configuration from the 
noninteracting valence state and its meaning can be clarified by considering the example 
of the Hydrogen atom. One can usually ignore the two-electron-one-positron configuration 
relative to the one with a single electron, not because two electrons plus one positron are 
much heavier than one single electron, but rather due to the fact that the two-electron-one- 
positron configuration is about 1.5 MeV (3m e ) off the energy-shell of the free state, whereas 
the off-shellness of single-electron configuration is only due to the binding energy, ~ 13 eV. 

In the most strict sense DR requires that the reduced theory be renormalizable in three- 
dimensions and that corrections to correlations calculated in this reduced theory be sup- 
pressed by powers of 1/T. 

It is well-known that only the first few terms of the reduced theory are renormalizable 
in three-dimensions, even when all degrees of freedom are bosonic |llj. In our particular 
case, we also find a three-dimensional renormalizable Lagrangian only up to one-loop order. 
Nonetheless, we can still derive an effective theory that is capable of describing screening 
physics with an accuracy better than the leading order. This kind of reduced theory in 
general contains higher dimensional operators and needs to be defined in some well-defined 
regularization scheme. The coefficients of the reduced Lagrangian are determined by requir- 
ing that, in the appropriate kinematic regime, the relevant one-particle irreducible graphs 
calculated in the reduced theory match the corresponding ones in the original theory. It 
is important to understand that, even though the coefficients are calculated perturbatively, 
the reduced theory is designed to maintain the infrared properties of the original theory and, 
hence, the solution of the reduced theory is in general nonperturbative [T3| . 



In the Gross-Neveu model |14| we have explicitly shown that the fermionic degrees of 
freedom that survive in the reduced Lagrangian are nonrelativistic. Futhermore, it has 
been suggested that the same happens to the quark degrees of freedom present in the 
QCD Lagrangian at high T. In this paper, we demonstrate that this is in fact the case 
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for QCD and hence derive the nonrelativistic Lagrangian for the quarks up to one loop 
with methods similar to those used in deriving the heavy-quark effective theory fl6|Jl7l . 
This effective Lagrangian should reproduce screening mass splittings up to order g A {T) and 
masses themselves up to order g 2 (T). 

We want to point out from the outset that our effective Lagrangian is relevant for de- 
scribing the long-distance behavior of spatial correlators only in those channels that are 
protected by some global symmetry (such as baryon number and flavor) from mixing with 
gluons, e.g., baryonic and mesonic screening correlators. Furthermore, this reduced theory 
is not meant for describing correlators that can mix with gluonic state (their large distance 
behavior is dominated by screening glueball masses), thermodynamical potentials or other 
bulk quantities that are dominated by the bosonic zero modes. In these later cases, the 
terms in the effective Lagrangian containing fermions give irrelevant contributions and can 
be dropped. 

The plan of the paper is the following. 

In the next section we outline the general strategy and the criteria for dimensional 
reduction when observables made with quarks are involved. Since the concepts we introduce 
in Sec. H are somewhat new and cannot be explicitly found in the existing literature, we make 
the discussion as complete as possible. In particular, we first review the screening physics 
and the observables the effective theory is intended do describe, e.g. spatial correlators, 
screening masses and other observables that are identically zero when quarks are explicitly 
ignored. The relevant kinematic region and the exact meaning of the mass-shell condition 
in the context of screening physics are specified next. We finally introduce the expansion 
parameter, the off-shellness, in the same section. 

In Sec. |l| we explicitly calculate the DR Lagrangian, both at the tree and one-loop 
levels in the specified kinematic region. Composite operators are considered in Sec. [TV]. 

In Sec. [V| we first discuss a power counting scheme that guarantees that the expansion 
parameter remains small and, therefore, that the expansion is self-consistent. Then we 
derive an effective Lagrangian for QCD, where quarks are treated nonrelativistically, in 
analogy with heavy quarkonium systems. The final form of the reduced Lagrangian and the 
relationships between parameters in QCD and in the reduced theory are explicitly given. 
We finally discuss the temperature regime where the reduced theory becomes quantitatively 
reliable. 

The final section contains the summary and our conclusions. 

II. GENERAL STRATEGY AND CRITERIA FOR DIMENSIONAL REDUCTION 

In this section we set the stage for generalizing the concept of DR to cases that explicitly 
involve fermions. First we give a brief review of screening physics and an overview of how 
dimensional reduction comes about: this intuitive guide to our calculation also serves the 
purpose of introducing an important concept, the "mass-shell" in screening quantities, that 
is extensively used throughout the paper. Then we state the criterion for DR to take place 
when observables whose leading contribution comes from fermions are present, after having 
recapitulated the corresponding criterion for DR in a pure bosonic case. Next we discuss in 
detail the relevant kinematic regions, relevant to screening physics, where DR could happen 
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both for fundamental fermionic fields and for composite operators. Finally, we outline the 
strategy for explicitly verifying whether DR occurs in QCD. Where appropriate we elucidate 
similarities and differences with the heavy quark expansion. 

A. Screening physics and mass-shell condition 

Finite temperature screening physics can be directly and naturally formulated in terms of 
a Lagrangian in the four-dimensional Euclidean space. Unlike real-time dynamics @, there 
is no need to eventually analytically continue results to the (3+l)-dimensional Minkowskian 
space-time, since screening phenomena are described by time-independent correlation func- 
tions. Taking advantage of this static nature it is convenient to Fourier transform the fields in 
the time direction in terms of Matsubara frequencies. Then the four-dimensional Euclidean 
Lagrangian can be equivalently rewritten as a three-dimensional Euclidean Lagrangian with 
an infinite number of Matsubara modes. 

We typically want to study correlations between operators in the limit of spatial distances 
much larger than the thermal wavelength, or more specifically, we only consider correlators 
at \x\ l/T. This large spatial separation selects a preferred direction, which we take 
along the first axis in the four- dimensional Euclidean space whose zeroth component is the 
imaginary time. Then the dominant large-distance contributions to such correlators come 
from the lowest singularities in the external momentum variable p\. 

If we consider a weakly interacting theory, it is intuitive that singularities appear when 
the external momentum is such that some of the denominators of the internal propagators 
vanish. Since the four-momenta are Euclidean and we consider massless particles, denom- 
inators are of the form M 2 + p 2 and vanish either when M 2 = and p 2 = (bosonic 
zero mode with vanishing spatial momentum) or when some of the spatial components 
is imaginary, p\ = —M 2 (bosonic non-zero modes M 2 = (2nirT) 2 and fermionic modes 
M 2 = ((2n — 1)ttT) 2 ). In this last case p\ is purely imaginary, which corresponds to an ex- 
ponentially decaying spatial correlation, a well-expected behavior when only non-zero modes 
are involved. 

In the complex plane of p\ the condition that the Euclidean four-momentum be zero 
becomes {ipi) 2 — p\ — p\ = M 2 , which can be interpreted as the mass-shell condition in a 
(2+l)-dimensional Minkowskian space-time for a particle with mass equal to the Matsub- 
ara frequency M. The concept of particles being on, close or off their "mass-shell" used 
throughout this paper to describe screening physics has the precise meaning given by this 
interpretation of the singularities in p\ in the original four-dimensional static correlators. 
It should be emphasized that the screening singularities we have introduced here have no 
direct relation to real-time dynamics, which corresponds to singularities in real frequency in 
the (3+l)-dimensional nonstatic correlators 0. 

Finally, it is helpful to remind here that thermodynamical quantities are determined from 
static Green's functions in the kinematic region with vanishing or small external momenta, 
in contrast with the screening physics directly involving quarks mentioned above. 

B. Criteria for dimensional reduction 
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1. Pure bosonic case 



DR in the bosonic sector can be explained by the observation that there exists a clear 
scale separation when temperature is high enough. Modes with non-vanishing Matsubara 
frequencies have masses of order T, while modes with zero Matsubara frequency have masses 
of order one. If one is only interested in the dynamics of the static modes in the low 
momentum regime, nonstatic modes might play negligible role, since their effects can only 
be felt through virtual processes with energies comparable to or higher than their masses. 

The above intuitive picture can be formalized within the framework of perturbation 
theory in the following way. The (D+l)-dimensional Lagrangian C D+1 is said to undergo 
DR to a specific .D-dimensional Lagrangian C D , in the high-T limit and to a given order in 
the coupling constant, if all the static Greens' functions of C D+1 , G D+1 (p,T), are equal to 
the corresponding Green's functions of C D , G D (p), up to corrections suppressed by powers 
of 1/T: 

G D+1 (p,m,T) = G D (p,m) + 0(\p\/T,m/T) , (1) 

where m represents possible external dimensionful parameters of the theory, such as a usual 
mass parameter. In general, the form and the parameters of C D are determined by the 
original theory. 

These naive expectations based on a tree-level power counting fail because of dynamically 



generated scales of order T |TT[ . Nevertheless, these dynamically generated scales only induce 
corrections proportional to powers of the coupling. If the coupling is small, the concept is 
still useful and the reduced theory still provides a simplified physics picture. 

It is apparent that if the original theory is renormalizable in D+l dimensions, the reduced 
theory is a super-renormalizable theory in general, since the coupling in the reduced theory 
will have positive dimension in mass units. However, whether the ultraviolet behavior of the 
reduced theory is truly improved depends on whether all Matsubara modes, except a finite 
number, decouple. Finally, we remind that beyond tree-level DR manifests explicitly in the 
reduced Lagrangian only in certain specific classes of subtraction schemes [§,0], although 
the physics of DR for a given theory is of course a scheme independent phenomenon. We 
shall make further comments on this point, when we discuss the choice of the coupling 
constant in the reduced theory. 



2. Fermionic case 

Due to the antiperiodic boundary condition in the temporal direction, all fermionic fields 
have tree-level masses of order T. As a consequence, fermions are usually treated like implicit 
degrees of freedom in the high-T limit. However, there are situations where we want to study 
observables that directly involve fermions and that are zero if no explicit fermionic dynamics 
is kept. A typical example is the correlation between mesonic currents, which is defined in 
terms of fermionic fields: the Lagrangian does not include terms that couple such currents 
directly to bosonic fields and the leading contribution is given by fermionic modes. When 
studying such cases, does it still make sense to ask whether some of the fermionic modes 



6 



are more important than others? The naive answer might be negative, since there exists no 
obvious scale separation between the lightest fermionic modes ±7rT and the rest. 

This apparent lack of a scale separation is also reflected in the fact that the typical 
spatial momenta of the fermionic modes at high temperature are not small relative to T. 
For example, even in the free theory case, fermions are always on their mass-shell and hence 
\p\ is of the order T. This implies that the expansion in \p\/T is meaningless. 

However, QCD is asymptotically free and quarks are weakly interacting in the high-T 
limit. As previously discussed, in this weakly interacting regime, quarks are almost on 
their "mass-shell". This intuitive picture suggests us that the relevant small scale is not 
the typical spatial momentum, which is large even in the free theory, but the amount by 
which the interaction brings the spatial momentum off the free-theory mass-shell (the "off- 
shellness" ) or, in other words, the residual momentum after the contribution from the mass 
has been subtracted out. 

If we consider, for example, one of the lowest Matsubara modes lo± = ±7rT, we can define 
the dynamical residual momentum q 2 = p 2 + (vrT) 2 and, similarly to the bosonic case, we 
say that the theory undergoes DR, if all the static Greens' functions G D+1 (p, T) are equal to 
the corresponding Green's functions G D (p) up to corrections suppressed by powers of 1/T: 

G D+1 (p, m, T) = G D (p, m) + 0(\q\/T, m/T) , (2) 

where again m represents possible external dimensionful parameters of the theory. Similarly 
to bosonic cases, corrections are suppressed only by powers of the coupling constants instead 
of powers of 1/T when there are dynamically generated masses proportional to T or, as we 
shall see, when the specific dynamics makes the residual momentum proportional to T. 

3. Analogy with the heavy-quark expansion 

The situation we have described for the high-T expansion in the quark sector shows 
several similarities with the heavy quark expansion in QCD |lq-|l8|. In fact, the heavy 
quark degrees of freedom can usually be integrated out leaving power suppressed corrections 
if our interest is in light quark observables, such as the D meson. When we also want to 
study observables that involve heavy quarks, such as heavy quarkonia, the explicit heavy 
degrees of freedom must be retained. In addition, the effective theory that describes the 
heavy quark sector is also derived by expanding in the residual momentum relative to the 
heavy quark mass. 

However, there are two major differences between the heavy quark expansion and the 
high-T expansion for quarks. One difference is that, in the latter case, we need to integrate 
out an infinite number of Matsubara modes for each flavor, whereas in the former only the 
antiquark degree of freedom is integrated out for each flavor. The other difference is that the 
coupling constant in the high-T reduced theory is, for pure dimensional reasons, proportional 
to T, while the coupling constant in the heavy-quark effective (HQE) theory, is independent 
of (or at worst logarithmically dependent on) the heavy quark mass. The first difference leads 
to the consequence that the HQE theory maintains the original dimensionality 3 + 1 and the 
high-T effective theory has its dimensionality reduced by one. The second difference implies 
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that, contrary to the HQE theory, the accuracy of the high-T effective theory is usually 
worsened from powers of 1/T to only powers of the coupling. 

C. Relevant kinematic region for fermion DR 

Following the strategy described at the beginning of this section, we examine here what 
are the precise kinematic regions where DR can take place in the fermionic sector for fun- 
damental and composite operators. 

From now on we call light modes the static gluon (Q n = 0) and the lightest quarks 
{uj n = ±7rT), while the rest of the modes are denoted as the heavy modes. 

1. Fundamental fields 

An external light quark close to its mass shell has a momentum of the form p = (cu±,p) 
with p 2 ~ — 7T 2 T 2 . For definiteness, let us consider the "particle" characterized by uo + : the 
same considerations can be trivially repeated for the other light mode u>-. Here and in the 
following each quark mode of different Matsubara frequency is regarded as a "flavor" in the 
three-dimensional theory. Quotes are used to distinguish this use of the word "flavor" from 
the usual one. Since the on-shell condition is defined in Minkowskian space, we choose pi 
as the "time" component or energy in the reduced world (the original time component p 
becomes a chirally invariant mass). Then the residual momentum is defined as q\ = pi — iu + , 
<?2 — Vi an d qz = Pz with \qi\ < T. The alternative choice of qi = pi + iuo + can be interpreted 
as the interchanging of particle and antiparticle of the same "flavor" . 

Of course, one can also choose pi ~ iu n T with n ^ ±1, i.e. pi being close to one of 
the heavy mass-shell. However, the physically most relevant singularities are those closest 
to the origin in the complex plane of p±, since only these are practically measurable on the 
lattice. Therefore, we only limit ourselves to the cases of p\ ~ iu±. 

A static gluon close to its mass-shell has a momentum of the form k = (0, k) with 
\k\ <C T. Therefore, the definition of the high-T expansion as an expansion around the 
mass-shell of the weakly interacting modes reproduces the usual condition when applied to 
the bosonic sector (pure Yang-Mills case). 

2. Composite operators 

The rationale underlying the choice of the relevant kinematic regime when composite 
operators, such as meson and baryon currents, are present is the same as the one used in 
the quark and gluon sector: we expect small deviations from the free theory. 

For a static mesonic current with momentum k = (0,k), the lowest singularity starts 

at k 2 (2ttT) 2 , as one can easily verify considering the free fermion-bubble graph. This 

singularity corresponds to the lowest particle-antiparticle threshold, where particles and 
antiparticles are defined as the excitations near the two possible mass-shell conditions: p\ = 
iuj + and pi = —iu + . To be consistent with this condition, the momenta of the particle 
and antiparticle in the composite operator have the form p = (u + ,p) and p' = (cu + ,p f ) 
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with pi ~ iu + and p[ ~ — iu + and k = p — p'. The momenta running along the free 
fermion-bubble show clearly why we call, in the reduced theory, quarks with the same po 
but "opposite" p 1 components particle and antiparticle, while we distinguish quarks with 
different p components using different "flavor" . 

We could perform a similar analysis for baryonic currents. Since we do not anticipate any 
further conceptual problem in this kind of extension, we consider in this work only mesonic 
currents. 

D. Expansion parameter and power counting 

1. Tree level 

In the relevant kinematic region, we can classify propagators according to their behavior 
close to mass-shell and then construct a corresponding power-counting scheme that charac- 
terizes the behavior of a given graph. For example, it is trivial to verify that, when close to 
mass-shell, a light-fermion propagator is of the order T/\q\, while a static gluon propagator 
is of the order of T 2 /k 2 , relative to propagators of heavy modes. 

The energy-momentum conservation at each vertex forces internal quark lines of a graph 
whose external lines are close to their mass-shell to remain themselves close to the mass- 
shell, unless at least one of the internal lines is heavy. However, a graph that involves 
internal heavy lines is suppressed relative to the same graph with the heavy lines replaced 
by light ones, as it is trivially demonstrated by the propagator classification given above. 
This result, in turn, implies that any tree graph with only light external quark or antiquark 
lines is correctly reproduced by the corresponding tree graph in the reduced theory. 

2. Loop effects 

The fact that the dominant contributions to screening observables come from the kine- 
matic regions where the external line are close to their mass-shell in the (2 +1) -dimensional 
Minkowskian space implies that the relevant power counting that establishes the relative 
importance of the difference graphs is not the usual one in Euclidean space. This necessity 
of explicitly considering the contributions of the Minkowskian singularities makes the power 
counting for graphs involving loops less straightforward than at the tree level. Fortunately, 
the number of graphs at a given order in the loop expansion is finite. At the one-loop level 
this number is small. We can just perform an explicit calculation and isolate the important 
contributions. 

The strategy we adopt is the following. We consider any graph, G4d{p, T), in the original 
theory in 3+1 dimensions, where p is close to the mass-shell, and we then decompose this 
graph into terms recognizable as three-dimensional graphs. We achieve this result by first 
separating out those terms in G iD (p,T) that involve only light lines: these terms are by 
definition reproduced by the reduced tree-level three-dimensional Lagrangian and we denote 
them by Gzd{p,T). The remaining contributions to G^oiPiT) involve at least one heavy 
line and we define them as G H (p,T). We say that DR occurs if all such G H (p, T)'s can be 
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made local, i.e. if they result in a polynomial in the residual momentum q. In other words, 
DR occurs, if 



G iD (p,T) = G 3D (p,T) + G H (p,T) = G 3D (p,T) + G H (0,T) + 0(\q\/T) , (3) 

where G (0, T) has either the form of terms already present in the tree- level reduced La- 
grangian (its effect is the renormalization of the relevant parameters) or the form of a new 
renormalizable vertex. More generally, the reduced theory needs also to contain nonrenor- 
malizable vertices if we want to reproduce graphs to higher orders; these terms can be treated 
consistently in the context of an effective theory (see Sec. |V|). 



III. THE REDUCED THEORY: CALCULATION 

A. The tree-level Lagrangian 

At the tree level the reduced theory can be simply obtained by Fourier transforming the 
QCD Lagrangian and then retaining only the static gluonic fields and the quark fields (after 
rescaling a factor of yT) with lowest Matsubara frequencies (u± = ±7rT), 

£rd = -\f%F% - \{D t A Y{D t A,y + + *MSy)7o + «7iA]tfo , (4) 

where F* = fyA? - djA" - g z f abc A\A), (DiA ) a = d { A% - g z f abc A\Al and D$ = - 
ig-^Af (A a /2)-?/>. For simplicity we have assumed that quarks are massless and have Nf flavors. 
Figure [I] shows the two vertices of £p R that involve quarks: the graphical notation is such 
that a thick (thin) solid line represents a quark with frequency uj + (o>_) and a wiggly line 
represents a static gluon. 

The coupling g 3 is related to the four- dimensional coupling through g\ = g 2 (fi)T. At the 
tree level the subtraction scale \i is not yet specified. Since £ RD is a super-renormalizable 
theory in three dimensions, all the dynamical scales must be set by the coupling constant 
g\ = g 2 (fi)T and temperature T. 

Of course, once loop corrections are included, the reduced theory in Eq. (P acquires new 
vertices and the coupling constant g\ has a more complicated dependence on the original 
coupling g 2 (/i). For example, g\ can receive corrections such as g 4 (/i)T and so on. However, 
since QCD is asymptotically free, the appropriately defined coupling constant (DR is only 
manifest in subtraction schemes that require \x ~ T [|TT1 . [19| , [14| ) has the asymptotic behavior 



g 2 {[i ~ T) ~ 1/ InT. Therefore, the corrections to the tree-level form of g 3 should not modify 
the fact that the two dynamical scales be g 2 ~ g 2 (fi ~ T)T and T itself at high temperature. 
The criterion according to which we choose \x ~ T is in general to minimize loop corrections 
to the leading result. We shall discuss the precise choice of the proportionality constant 
later in the paper. 

The vertices that involve light quarks in the reduced theory, the ones shown in Fig. [I], 
are not the only vertices present in the original theory. There are additional vertices that 
involve at least one heavy mode; we show them in Fig. [2] with the notation that a double 
line stands for a heavy quark mode and a spring-like line stands for a nonstatic gluon mode. 
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These vertices are collectively called A£h- In the first vertex the fermionic line does not 
change its "flavor" , since the gluon is static, and we call it "flavor-conserving" . All the other 
vertices in Fig. involve a nonstatic gluon and, therefore, the fermionic line changes its 
"flavor" ; we call these vertices "flavor-changing" . 

Since we are only interested in graphs with light modes in the external lines, energy- 
momentum conservation implies that these heavy vertices cannot contribute at the tree 
level. Our job is to verify whether the corrections induced by these heavy vertices at the 
one- loop order can be accounted for, in the high-T limit, either by readjusting the parameters 
of £q R or by adding additional local and renormalizable vertices. 



B. One-loop graphs 

In the following calculations we use the dimensional regularization for spatial integrals, 

r Jih, oo „ j3-2ejU oo „ 

and the MS subtraction scheme. For convenience, we work in the Feynman gauge and use 
the Euclidean Feynman rules given by Ramond p0| . 



The subtraction point fx needs to be proportional to T to make the DR manifest ||Tl 
then we define the T-dependent coupling constant 

g\T) = g* m (fx) (6) 

by choosing 

fi = (Aire-^-^T . (7) 

The temperature independent constant c can be read as a convenient way of parameterizing 
one of the possible subtraction schemes. The specific choice c = (N/2 — 2iVf In4)/(11./V — 
2Nf) corresponds, for instance, to the scheme that makes DR optimal for the pure-gluon 
effective action in the background field method We shall comment later on alternative 



choices such as the one that makes DR optimal in terms of the quark-gluon vertex function. 
In the following we also use the auxiliary coupling 

g\T) = ~g 2 (T)/(16n 2 ) , (8) 

for the sake of making formulae more compact. 

Without any loss of generality, we only show results for cases with external quark-line 
frequency u = u + = ttT; trivial modifications yield the corresponding results for the other 
light "flavor" with u = uo_ = —ttT. Moreover, we select the particle sector, which we have 
defined as the kinematic region where the first component of the spatial momentum is close 
to again there are only trivial differences for the antiparticle sector, which is defined by 
the alternative choice of the spatial component being close to —iuj + , and that is separated 
from the particle sector by a off- (mass) shellness 2ttT. 
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In the Matsubara-frequency loop-sum, the term with zero frequency can be easily rec- 
ognized as the contribution from the three-dimensional tree-level Lagrangian itself, up to 
trivial factors of T. Therefore, when computing the loop corrections due to the heavy Mat- 
subara modes, we leave out the term with n = 0, which is the direct contribution from the 
light modes and it is already generated by the reduced theory. 

The calculation of the one-loop amplitudes and their high temperature expansion around 
the appropriate kinematic regions can be done following a standard procedure, which, in 
general, involves the following steps: 

(1) combine the denominators by means of the Feynman parameter representation and 
perform the spatial momentum integral; 

(2) shift the external momenta: p = q + iuj + e\ and p' = q' + iu + ei, 

(3) expand the result in terms of the residual momenta over T, e.g. \q\/T; 

(4) perform the integrals over the Feynman parameters; 

(5) perform the Matsubara sum and express the result in terms of the Riemann zeta function. 
In the following we only give the final results obtained by the application of the above- 
described procedure. 



1. Quark self-energy 

The quark self-energy correction due to heavy modes is found by calculating the graphs 
drawn in Figs. [3] (b) and (c) with external momentum p = (u + ,p) = {u +) iuj + + qi, q 2 , (fe). 
The result is: 



is(p) = -c f g 2 (T) | [(x + 2c) p • 7 + x^ lv ] + I 



(9) 



where Cf — (N 2 — 1)/ (2N) and the coefficients X and X^ v are pure numbers that are listed 
in Table |. In this Table there appear the derivative of the Riemann zeta function evaluated 
at —1, £'(— 1), and the Euler's constant gamma je, whose approximate numerical values are 
C'(-l) « -0.16542 and ^ E ~ 0.57722. 

We find that the heavy mode correction to the self-energy E(p) is suppressed relative to 
the tree- level piece ip-j by the factor Q 2 {T). In addition, we also point out that, as expected, 
no chiral-symmetry-breaking mass-term for the quark self-energy has been perturbatively 
generated, in spite of the fact that other noncovariant terms have instead appeared. 



2. Quark- gluon vertex 

The corrections of the heavy modes to the quark-gluon vertex come from two types of 
graphs: graphs that have an analogue in QED, Figs. f| (b) and (c), and graphs that are 
intrinsically nonabelian, Figs. § (e) and (f). In the following the momentum of the incoming 
quark is labeled by p = (u + ,p) = {oj +1 iuj + + qi, q 2 , (fe) and the momentum of the outgoing 
quark by p' = = (u+, iuj+ + q' 1 ,q' 2 , <&). 

The first type of graphs yields: 
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r,(p,p / ) = i ~g(T)-(C f 




(Y + 2c) llM + Y^ lu ]+0( 




(10) 



where C a d = N and the coefficients Y and Y^ u are given in Table |. 
The second type of graphs yields: 

r>,p') = i~9(T)^rC ad G 2 (T) I \(Z + 3c) 7^ + Z^ lv 




(11) 



and again the coefficients Z and Z^ are listed in Table |. 

In Table | one can notice that there exist precise relations between some of the entries 
for the quark self-energy (X's) and of the abelian part of the quark-gluon vertex (V's): the 
reason of these relations is the existence of generalized QED-like Ward identities at finite 
temperature, due to the static gauge invariance. The fact that our results verify these Ward 
identities serves as a very useful consistency check. 

We find again that the heavy mode corrections to the tree-level coupling coming from 
T[i(p,p') and T' {p, p') are suppressed relative to the tree-level piece —ig{T)^^ by the factor 
Q 2 (T). In addition, there are also new vertices not present at T = 0, but these vertices are 
also order of Q 2 (T) relative to —ig(T) r ) tl . 



3. Vacuum polarization 

The contributions to the vacuum polarization tensor coming from both light and heavy 
quarks are given by the graphs shown in Figs. |5] (a), (b) and (c). There are a few small 
differences relative to the previous two cases. One is that we can consider the three graphs 
together and we do not need to separate out the contribution from the lowest modes, since 
all three are infrared finite. In addition, the external gluon line is static k = (0, k) and the 



loop Matsubara sum is now fermionic p 



[UJr 



, p), with u n = (2n — 1)itT. 



In the magnetic sector ITjj(fc) remains transverse 

IL,(fc) = (kikj - k%)8 ab N f 9 ^{l + lE + ln(i^) + O(^)} . (12) 

After renormalizing at the scale \x = 47re~ 7 ~ c T, we see that the only effect of is to give a 
finite wave-function renormalization to the static gluon compared to the zero temperature 
case. 

We also find that Ho = Hoi = 0, which in turn implies that there is no mixing between 
electric and magnetic components in the static sector induced by the nonstatic modes. 
The result for the electric part is 



rioo(fc) 



SabNf 



9_M 

24tt 2 



4tt 2 T 2 + k' 



- +7£ + hi(-^ - 1 



(13) 



After the wave-function renormalization, which again only receives a finite contribution with 
respect to the zero temperature case, we are still left with an additional T 2 term that implies 
a mass generation for the time component of the static gluon field, which is the well-known 
analogue of the Debye screening in QED. 
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C. Summary of one-particle irreducible graphs 

The effect of the heavy modes of the original four-dimensional Lagrangian, which are 
no longer present in the reduced theory, can be reproduced, at the one-loop level and in 
the relevant kinematic region, by the following three types of corrections to the tree-level 
(2+l)-dimensional reduced Lagrangian. 

(a) The quark bubble insertion corrects the static gluon propagator and is suppressed by a 
factor Q 2 {T). It also generate a mass term of order Q 2 {T)T 2 for Aq. 

(b) The quark self-energy insertion corrects the lightest quark or antiquark propagator and, 
in addition, generates new terms that are not present at zero temperature. Corrections and 
new terms are all suppressed by a factor Q 2 {T). 

(c) The static gluon-quark vertex insertion corrects the tree-level vertex and generate new 
vertices. Again corrections and additional vertices are suppressed by a factor Q 2 (T). 

(d) In this section we only considered graphs explicitly involving quark fields. Since the 
on-shell condition for bosonic fields are the same with or without the presence of fermions, 
graphs containing purely gluons or ghosts have already been considered in Refs. [^T],|T2J. 

In summary, we find two kinds of corrections in the infrared limit, namely when the 
light modes are close to their mass-shell. There are corrections that are directly generated 
by the corresponding one-loop graphs in the reduced theory, once the coupling constant 
is properly chosen. In addition, there are also one-loop corrections that are not contained 
in the tree-level reduced theory. However, all these terms are infrared finite to the order 
considered and are suppressed by a factor Q 2 {T) and hence subleading. The fact that this 
corrections are infrared finite implies that they can be accounted for by adding new local 
and renormalizable vertices to the reduced theory. 

IV. COMPOSITE OPERATORS 

In confined theories, such as QCD, the fundamental degrees of freedom are not manifest 
in the spectrum and it is necessary to use composite operators to probe physical particles 
of the theory. For example, one uses composite operators as interpolating fields for mesonic 
and baryonic states. Therefore, the study of how composite operators are reproduced in the 
reduced theory is necessary to have a complete picture of the DR physics. We shall find 
that considering composite operators introduces new features that are not trivial extensions 
of what already discussed. 

For the sake of concreteness, we focus our attention on flavor nonsinglet mesonic currents. 
Generalization to other cases, such as baryonic operators, can be done analogously. At the 
tree level, the static limit of these currents can be written as a sum over Matsubara modes 

O v = fdr i>(r, x)T^(r, x) oc ]T ^ n (x)T^ n (x) , (14) 

J n 

where T is any of the sixteen Dirac matrices. As discussed in Sec. O, the kinematic region 
of interest is the one where the lightest modes are close to their mass-shell: 2ttT in mesonic 
and 3nT in baryonic cases, respectively. Then at the tree level the high temperature limit 
implies that the dominant contribution comes from the operator obtained by using only the 
lowest Matsubara quark modes. 
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A. One-loop correction 



The procedure for calculating the one-loop correction to composite operators is similar 
to the one used for calculating the one-loop vertex correction, except that the momentum 
carried by the composite operator is now close to the particle-antiparticle mass-shell, in 
contrast with the momentum carried by a static gluon in the vertex correction. The explicit 
graphs are shown in Figs. ^ (a), (b) and (c). The choice of the flavor nonsinglet current 
avoids the mixing with gluonic fields^. 

More specifically, the external quark and antiquark momenta are expanded according to 
p = q + and p' = q' — iu + ei, where q and q' are supposed to be small relative to T. 

This choice reproduces the expected on-shell condition for mesonic currents p — p'm 2iuj + e\. 
Since this kinematic difference implies a new feature, we give more details in this case. 

The one-loop correction to the composite operator is proportional to 

AO r = 9 2 (v) C f V^ v {p,p) 7a7/J7^7a , (15) 
where V^ y {p,p') is defined as 

V, u (p,p) = 2T £ / [dk] I" da C' a d(3 (P + kW + k) v 

n& J Jo Jo (l-a-p)k 2 + a(p + k) 2 + {3(p> + k) 2 



This definition obviously implies that V^ u (p,p') = V u ^{p',p). The term in the Matsubara 
sum with n = 0, which is represented in Fig. || (a), has been excluded from V^ u , since this 
term is directly reproduced by the reduced theory £^ D , given by Eq. @. 

Fig. ^| (c), i.e. the terms with n ^ —1,0, can be checked to be infrared finite; their 
explicit contribution to V^ v (p,p') is, up to corrections of 0(\q\/T, \ q\/T), 



Vr 



M = ^{7 + H^) + ™ + w ™} . ( 17a ) 

^ = + (17b) 

4 C) (p, P) = "4 C) (P, P) = - W01 , (17c) 
Vi c) (p, p') = -^ c) (p, p') = for i = 2, 3 , (17d) 

where the numerical coefficients W 00 , Wn, and W i are given in Table |. After taking 
care of the 1/e terms, which have exactly the same form as at T = 0, using the standard 
composite operator renormalization, there only remain the infrared finite terms that can be 
reproduced by correcting the tree-level current Or in the reduced theory with the addition 
of local operators. 



1 States that mix with gluonic fields are dominated by the bosonic zero modes and can be described 
by the pure glue reduced theory with the addition of local operators that couple quarks and gluonic 
states. 
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B. Additional infrared singularity 



Now let us focus on the remaining term shown in Fig. [| (b), corresponding to n = — 1 
in Vfj, u (p,p') of Eq. fll6|) . In this term a heavy gluon of frequency Q n = —2nT is exchanged 
with a corresponding change of the "flavor" (frequency) of the quarks at the vertices. The 
expansion around the mass-shell of this term yields, apart from terms of 0(\q\/T, \ q'\/T), 

^W) = y^x (18) 
j [ X ~°ja + ^oM 1 -a + p)- z^iM 1 +a-/3) + 5^5 vl [l - (a- P) 2 ] 



l rL—a 

da / dp 
Jo 



o A) [(2 - a - P) 2 + 2aa(l - a) + 26/3(1 - p) - 2ap(2 -a-b)} 



3/2 



where a = iqi/(nT) and b = —iq' 1 /(irT). The Feynman-parameter integrals can be carried 
out with the result 

+ 32^ + ^ 41n2 - 3 )] + ^) ' (19) 

which is clearly logarithmically divergent in the infrared. 
The situation is summarized as follows. 

(1) The term represented by the graph in Fig. ^| (a) diverges linearly and gives the leading 
infrared contribution coming from the composite operator: this term only involves the lowest 
Matsubara modes and, therefore, its infrared physics is exactly reproduced by the reduced 
theory C^ R . 

(2) The terms represented by the graph in Fig. |] (c) are infrared finite: these terms involve 
heavy Matsubara modes and their contribution is not generated by the tree-level reduced 
Lagrangian £r D , but they can be compensated by adding new local operators to Or- 

(3) The term represented by the graph in Fig. [] (b) diverges logarithmically and gives 
a subleading infrared contribution compared to graph (a): this term involves the lowest 
Matsubara modes for the quark lines and a heavy mode for the gluon line and, therefore, 
it is not present in the reduced theory. In addition, this infrared logarithmic behavior 
apparently implies that it cannot be generated by adding a local correction to the operator. 
In other words, we need to settle the question whether it is possible to add to the tree- 
level operator in the reduced theory higher- dimensional operators that might generate the 
same logarithmic singularity. This option is easily ruled out, once one recognizes that 
this logarithmic singularity is associated with the total external momentum carried by the 
composite current q% — q[ and is related to the particle-antiparticle (of the "— " mode) 
production threshold, as can be seen explicitly from Fig. |] (b). 

We are left with the option of adding a higher-dimensional vertex to the reduced theory. 
We can easily guess the form of this new vertex from the fact that the logarithmic singularity 
comes from the kinematic region of Fig. || (b) where the particle and antiparticle are close to 
their mass-shell, while the heavy gluon is far off its mass-shell, i.e. the gluon line "contracts" 
to a point. In fact, integrating out the heavy gluon yields the four-quark vertex 

AC F = f^ + 7,yV>-) , (20) 
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which is depicted in Fig. ^ (a). If we add this four-quark vertex to the reduced Lagrangian, 
there appears a new contribution, shown in Fig. [7] (b), to the composite operator one- loop 
calculation. This new contribution gives the same logarithmic singularity of Eq. (|T9| ) 



ln(a + b) 
64tt 2 



+ 



32tt 2 



(^o^o + i$tio$vi - iS^iS u0 + S^ui) In 2 - {8^8 jX + 8 i3 )j2 



(21) 



while the finite differences can now be compensated by local corrections to the operator. 

Of course, if we add this new four-quark vertex to the reduced Lagrangian, we need also to 
consider its additional contribution to the fundamental one-particle irreducible graphs which 
have already been calculated in Sec. [TIT]. In particular, we must check that the corrections 
to those graphs induced by this four-quark term are infrared finite. At the one-loop level, 
we need to consider only two graphs: the one depicted in Fig. [7| (c), which contributes to 
the fundamental vertex, and the one depicted in Fig. |7] (d), which contributes to the quark 
self-energy. An explicit calculation shows that these two contributions are indeed infrared 
finite. Explicitly, the correction to the vertex, Fig. |7| 



is 



A 6 



T,(p,p') = -ig(T)-(C 



a 



ad 



G 2 {T) <7o7m7o 




(22) 



while the correction to the self-energy, Fig. [7| (d), contributes to the chirally invariant mass: 



it{p) = C f g\T) 2uo +l0 



(23) 



Finally, one might worry that the four-quark term given in Eq. (P0f), since its mass 
dimension is four, might be nonrenormalizable in the reduced theory, at least by a naive 
power-counting argument. However, the use of dimensional regularization makes all graphs 
shown in Figs. [?| ultraviolet finite. In fact, it is the very fact that this operator has mass 
dimension four that makes its contribution, which is suppressed by one power of 1/T at the 
tree level (AC F ~ 1/T), be actually only suppressed by g 2 (T) in the one-loop graphs of 
Figs. (b), (c) and (d). 



C. Comments to the one-loop calculation 

In principle, one could go on and examine higher orders in the loop expansion. However, 
this additional effort is useless in this framework. In fact, the naive expansion in powers 
of 1/T is only possible, in general, to the leading order in g 2 (T). First, thermal masses 
are eventually generated at some order in the coupling and these masses break the implicit 
assumption that all the relevant momenta can be made small at will. Second, the power 
counting argument that is discussed in the next section shows that the residual momenta 
(off-shellness) in quantities such as screening masses is of the order q 2 ~ g 2 (T)T 2 , which 
implies that terms proportional to q 2 /T 2 are not suppressed by power of 1/T but only by 
powers of g 2 {T), which vanishes only logarithmically. 



17 



At this point we could write down the reduced theory at the one-loop level, utilizing the 
calculations done earlier in the last and this sections. However, one important point about 
the expansion of the off-shellness is the following. When we choose the first component of 
the momentum (pi) to be close to either ±iirT as the starting point of the off mass-shell 
expansion, this choice breaks the symmetry between what we call particle and antiparticle 
in the reduced theory. If we insist that the reduced theory be formally relativistic, this 
asymmetry results in the entanglement of corrections of different orders in g 2 (T) at a given 
order in loop expansion. In other words, the heavy quarks of the (2+l)-dimensional reduced 
theory become nonrelativistic (NR) in the high-T limit, as we shall see in the next section, 
and the relativistic formalism has the effect of retaining degrees of freedom that are higher 
order than the ones already dropped. Therefore, the achievement of an accuracy better than 
the leading order requires the systematic separation of the particle and antiparticle contri- 
butions, together with the correct counting of the contributions from the "flavor- changing" 
term of Eq. (POD- We then postpone the explicit expression of the reduced Lagrangian to 
the next section, where we discuss this more systematic nonrelativistic reduction. 



V. NONRELATIVISTIC EFFECTIVE THEORY 

The possibility of describing the high-T QCD screening physics by a renormalizable local 
Lagrangian in 2+1 dimension, is valid only up to one-loop level. However, since our specific 
goal is to study the screening physics, it still makes sense to give up the renormalizability and 
derive an effective action for solving screening states in the high-T limit with accuracy better 
than the leading order. In fact, the effective action necessarily contains higher dimensional 
operators. The coefficients of these higher dimensional operators will be derived by using 
the so-called matching technique, in close analogy with the nonrelativistic reduction applied 



to the positronium |L6[ in QED and to heavy-quark systems |17| in QCD. 

In this section we derive a (2+l)-dimensional effective theory that describes screening 
states with accuracy up to g A (T), apart from an overall additive zero-point energy which we 
determine only up to g 2 {T). As discussed at the end of the previous section, this expansion 
requires an explicit separation of the particle and antiparticle sectors, because of the asym- 
metry of the mass-shell condition. Therefore, this new theory, which improves results for 
screening mass splittings by one order in g 2 (T), becomes necessarily nonrelativistic-like for 
quarks. 

The basic strategy is the following. We write down the most general nonrelativistic 
Lagrangian, with terms up to some power in the appropriate power counting scheme. Then 
the coefficients of these terms are chosen so that they reproduce Green's functions of the 
original theory expanded up to the same power in the coupling in the relevant kinematic 
region. 



A. Notations 

Since we are interested in the regime of the reduced theory where quarks are close to 
their mass-shell, in the sense discussed in Sec. ||, it is convenient to explicitly rotate from 
an Euclidean notation (lower case letters) to a Minkowskian notation (upper case letters). 
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We rotate the original first spatial direction to the time direction of the (2+l)-dimensional 
theory, while the other two spatial directions remain the spatial directions of the (2+1)- 
dimensional theory: from now on, bold face letters indicate these two-dimensional vectors. 
We label the original time direction as the third axis, which now represents the chirally 
invariant mass in the reduced theory. Specifically, the momentum of a quark mode with 
tree-level mass M = ttT is rotated according to 

(po = iM, Pl ,p 2 ,p 3 ) -> ((M + P ) = -i Pl ,Pi = pa, P 2 =p 3 ,M = p ) = ((AT + P ), P, M) . 
The related rotation of the Dirac matrices is 

(70, 7i> 72, 73) (r = iii, r x = -72, r 2 = -73, r 3 = -70) = (r , r, r 3 ). 

We use for the matrices Y the explicit form of Itzykson and Zuber [21|| : Yq — u 3 ® I 2 x2i 
r = ia 2 (g> cr, and T 3 = i<7 2 <8> ct3- This choice of the Y matrices yields the following explicit 
form of the on-shell spinors, which are defined by [E(P)Yq — P ■ T — MY 3 ]U(P) = and 
[E(P)Y - P • r + MY 3 )V(P) = 0, 



U(P)= 1 (<r-P + **M\ 
U[ ] V2E(P) { E(P) ) 




v(p) - _ :r. „ , (24) 



where E(P) = ^/M 2 + P 2 is the on-shell energy of a free quark in the reduced theory. 
These spinors specify the basis in which the relativistic 4-component quark field is decom- 
posed as the nonrelativistic 2-component (spin up and down) particle field and 2-component 
antiparticle field. 

For convenience, we also rename the original gauge fields in the reduced theory as: 
A\ = Aq, (A2, A 3 ) = A. and A = ip. Since in the reduced theory ip transforms like a matter 
field under the original static gauge transformation, it is often called the "Higgs" field. 



B. Dimensional analysis 

In general, the expansions in powers of g 2 {T) and in derivatives (or low momenta) are 
two independent expansions. However, when we consider bound states (screening states in 
our case) that are dominated by the perturbative interaction, the two expansions become 
intertwined. The reason is that the typical momentum is no longer an independent variable, 
but rather it is determined by the interaction, in contrast with scattering experiments where 
one controls the momenta externally. If the interaction is weak, the typical momentum is 
proportional to some power of g 2 {T). 

It is possible to develop a systematic method for counting the contribution of each term 
in powers of the coupling constant to this combined expansion. The rationale of this method 
is described in detail in Ref. |17ll , and it is basically based on the analysis of the relevant 
Schrodinger equation with a potential derived from the tree-level approximation: in our case 
V(x) ~ g 2 T\n\x\. The resulting power-counting rules valid for studying screening states 
(bound states of the reduced theory in 2+1 dimensions) are shown in Table 0. 

It is important to emphasize that the power counting is determined by the leading be- 
havior of the potential in terms of the coupling constant. Therefore, it is still valid even 
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when the potential acquires a nonperturbative linear confining term. In fact, in the high- 
T limit the perturbative tree-level potential overwhelms the induced spatial string tension 
e s (T)<x~g\T)T 2 @. 

Quarks in the reduced theory are very heavy making the time direction special relative 
to the spatial ones: this fact is reflected in the form of the effective Lagrangian, which is 
nonrelativistic, and also in the power counting rules of Table [II 



C. Tree level 

At the tree level there are two kinds of corrections to the tree-level NR Lagrangian: 
kinematics corrections and corrections to the elastic scattering of a quark from external 
sources. 

We do not need to consider inelastic scattering, because we impose that the reduced 
theory contain only one power of the time derivative d t : this requirement corresponds to 
the precise choice of the field parameterization given in the U and V basis. There are no 
inelastic terms at the tree level in the original theory, and the inelasticity only appears at 
the one-loop level with this specific choice. In principle, one can relate this particular choice 
of the field parameterization to others that involve higher power of dt using the invariance 
of the physics under field redefinition. 

According to the power counting rules shown in Table [H], the only corrections up to order 
g 4 (T) relative to the leading term are the following ones, 
fa) Kinematics correction: 



U(P)P ■ T U(P) — M + P - E(P) = P 



p2 p4 

+ 



2M 8M 3 



(b) Scattering from external Aq\ 

U(P)A (P - P')ToU(P') 
1 



Ao(P-P') 



AM 2 

[c) Scattering from external .4.: 



1 -(P - P'Y + -a 3 P x P' 



2M 2 



2M 



<x x (P - P') 



U(P)A(P — P') T U(P') = Ai{P - P') 
(d) Scattering from external (p: 
U{P)(p{P - P')T 3 U(P') = y3(P - P') 



— (P l + Pi) - — e i7 (P,- - P> 3 
2M V lJ 2M n J ] ' 



_L(p2 + p 2 ) + — ^ — (j x (P - P ) 

4M 2 K J 2M V 1 



(25) 



(26) 



(27) 



. (28) 



D. One- loop level 

In Sec. |II| we have already calculated the heavy mode contributions at one loop. Since 
now we are using a NR formulation of the reduced theory in the quark sector, the light 
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quark contributions of the original theory are not exactly equal to the ones in the NR 
reduced theory, even if they have the same infrared behavior. Therefore, we also need to 
calculate the one-loop contribution of the light quark modes both in the original theory 
and in the NR reduced one: the difference between the two results need to be added as 
a correction to the reduced Lagrangian together with the heavy mode contributions. In 
addition, the one-loop amplitudes needs to be sandwiched with the appropriate spinors to 
yield the correct correction terms. The calculation is straightforward and we only list the 
resulting terms up to p 4 (T). 

The nonrelativistic one-loop corrections, which need to be compared to the analogous 
corrections of the light quarks in the original theory, are calculated using the following 
tree-level Lagrangian 



NR 
tree 



ifr 



iD t + 



El 

2M 



if) 



(29) 



where D t = d t + igzAo, D = d — ig$A and if) now is a two-component (representing spin up 
and down) nonrelativistic quark field. One can easily work out the Feynman rules associated 
with this Lagrangian and calculate the relevant one-loop graphs. Again, we only list final 
results, which have been computed using dimensional regularization. 

In the following the one-loop corrections to the original theory both from heavy and light 
quarks are grouped together (the ones coming from sandwiching with appropriate spinors, 
the terms calculated in Sec. |TTJ and the "new" ones coming from the lowest Matsubara 
frequency). The one-loop contributions from the NR reduced theory are instead given sep- 
arately. The explicit values of the X's, Y's and Z*s are given in Table |T[ 
(e) Self-energy corrections: 

U(P)iZ U(P) = C f g 2 (T)\ -4(M - P ) 



+ 



2M - (X + X n - iX 10 + 2c) P + (X - X 00 - zX 10 + 2c) 



2M 



iZ NR = C f g' 2 (T)i-2M 



1 , 71 Ll 



3P^ 
4 ~M 



(30) 
(31) 



The momentum independent divergence in iTj NR is related to the additive mass renormal- 
ization in the nonrelativistic Lagrangian. This additive mass term is often explicitly ignored 
at the price of introducing a zero-point energy ambiguity in the nonrelativistic theory, i.e. 
we can only calculate mass differences. In principle this ambiguity can be resolved at least 
perturbatively in each specific regularization scheme, 
(f) Abelian-like vertex corrections: 



U(0)Ar U(0)=i~g(T)^(C f C 



ail 



2 

U(0)AT 3 U(0)=zg(T)^-(C 



2 

C a d 



-4 



Y + Y u - iY 



10 



at 



NR. 



AT 



NR 



g\T) 

)g\T){-A+ Y + Y 00 + tY 01 

2g 2 (T) . 



2c 



2c 



(32) 
(33) 
(34) 
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(g) Nonabelian-like vertex corrections: 



C/(O)Ar £7(0) = ig(T)¥- C ad Q\T) 1-2 + [z + Z u - iZ 10 + 3c 



£7(0) Ar 3 £7(0) = z$(T)y C ad Q\T) 



Z + Z 00 + %Zm + 3c 



AL^ = Ar^ K = 



(35) 

(36) 
(37) 



Since the leading terms are already order g 3 (T), the momentum dependent terms are higher 
orders and can be dropped in vertex corrections. 

As expected, the one-loop corrections in the original theory and the NR reduced theory 
do not match exactly. The differences can be compensated up to order g 4 (T) by adding the 
following terms to the (2+l)-dimensional NR Lagrangian 



f d 2 

A£i-i oop = Q 2 (T) | ci ipUdtip + c 2 ^VjjrV' ~ c 3 i> 9s AV' - c 4 W93<pip 



(38) 



The coefficients Cj's, which are determined by subtracting the one-loop corrections obtained 



in the NR reduced theory from those computed in the original theory, are given in Table |T] . 

It is important to realize that we have taken out a factor of 2 from the one-loop results 
computed in the original theory before we subtract the nonrelativistic results from them. 
The necessity of dividing out this factor of 2 comes from the fact that the Lagrangian in 
Eq. (HD has two degenerate "flavors" (with M = ±7rT) , whereas the NR one-loop results are 
for a single "flavor". If the physics requires a normal flavor structure, the relevant flavor 
content must be added explicitly to the NR Lagrangian. 



E. Reduced Lagrangian: fermionic part 

We can now write down the final result for the fermionic part of the DR Lagrangian by 
collecting the corrections from the last two subsections and properly gauging each derivative. 

The final form of this Lagrangian is more suggestive if written in terms of the appropriate 
color electric and magnetic fields (i = 1, 2) 

£i = Foi , £>3 = -^ij3^ij , (39) 

where T%j = diAj — djAi — igs[Ai, Aj] is the gauge field strength tensor in 2+1 dimensions. 
Note that there is only one component of the color magnetic field in 2+1 dimensions. The 
final Lagrangian is given by 

C F = £<°> + + + , (40) 

with 

£ (0) = (i Dt + ^) * - V *V* , (41a) 

£(1S) = "2%* to " X £ * + 2M* t<T X ( D(p )* (41b) 
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£(2) 
C (2s) 



93 

AM 2 



9s 
2M 



£ ■ DW 



AM 2 



D 4 



1 ' y/ 8M 3 



(41c) 
(41d) 



where the sum over the color indices is implicit. In the above equation we have absorbed the 
one-loop corrections, due to A£i_i oop in Eq. (|58D, into finite renormalization of the proper 
physical quantities: 



M 

93 

9 V 



1 + |^ 2 (T)]^ 



TTT 



l + ( Cl -c 2 )£ 2 (T) 
~g(T)Vf[l - ( Cl - c 3 )g 2 (T) 
~g(T)VT\l-( Cl -c,)g 2 (T) 



(42a) 
(42b) 
(42c) 
(42d) 



At this point it is appropriate to make some general remarks concerning the fermionic 
effective Lagrangian. 

(1) The Lagrangian in Eqs. fl4Tf ) and fl4T|) has contribution only from a "single-flavor" particle 
state. The reduced Lagrangian for the antiparticle state has the same form except that 
couplings change sign. The four-fermion term in Eq. (p0|) does not give contribution to 
order g 4 (T) because of the separation of the quark and antiquark sectors. 

(2) According to the power counting rule in Table [TIL begins to contribute to the binding 
energy at order g 2 {T)T } at order g 3 (T)T, and £ (2) and £ (2s) at order # 4 (T)T. While 

and are spin independent, and are spin-dependent. 

(3) Even though there exist NR four-fermion terms that are order g 4 (T) according to a naive 
application of the power-counting rules of Table O, e.g. (7 2 (^^ $"!"$) jM ($ is the antiquark 
field), these terms can only contribute to the binding energy through higher- loop graphs 
and, therefore, their contributions are in fact of order higher than g A {T\ 

(4) Notice that is absent in 3+1 dimensions since it is not a scalar under spatial 
rotation, while it is a scalar under a two-dimensional rotation around the 3rd direction. 

(5) Since the self-energy in the NR theory is specified up to an additive constant, the 
Lagrangian in Eq. (f40|) cannot give the zero-point energy. However, it is still possible to 
determine perturbatively the zero-point energy shift between the NR Lagrangian and the 
original one. For instance, the formulae given in Eqs. (|30|) and (|31|) give the zero-point 
energy shift up to order g 2 (T). If we want to reach the same accuracy we have obtained for 
the energy differences (splittings), i.e. g 4 {T), also for the zero-point energy shift we need to 
perform a two-loop calculation of the quark self-energy. 

(6) The fact that g 2 (T) is small, as we shall see shortly, means that corrections from heavy 
modes and, therefore, the coefficients of higher dimensional operators are small and can be 
calculated perturbatively. It does not mean that the physics governed by the DR Lagrangian, 
whose coupling constant g 2 (T)T is very large at high T, is perturbative. In fact, the infrared 
behavior of the DR Lagrangian, by construction, remains the same as the original theory. 

(7) The coefficients in Eq. (|40|) has been derived in the MS scheme. If one wants to solve 
the reduced theory in a different renormalization scheme, e.g. on the lattice, one needs, in 
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principle, to compute again these coefficients in that specific scheme. In practice, this may 
not be necessary, since Q 2 {T) has turned out to be numerically very small and, hence, the 
tree-level coefficients, which are scheme independent, dominate. One could have problems 
only in those scheme that have corrections anomalously large. 

(8) The full gauge invariance in the original theory is now reduced to the "static" gauge 
invariance, which is explicitly kept in the above Lagrangian. 

(9) In analogy to universality in the pure bosonic case, a similar universality also holds in the 
quark sector. After the one-particle irreducible graphs have been matched in the relevant 
kinematic region and the reduced Lagrangian obtained, any other Green's function can be 
computed from this same Lagrangian to within the same accuracy of the matching. For 
example, we expect that the reduced theory is capable of describing the screening states, 
e.g. those states that describe the spatial correlation of mesonic currents; these states can 
be interpreted as bound states in the reduced world and hence they are beyond those 1PI 
graphs we have used to derive the reduced theory. We have explicitly verified this fact in 
the Gross- Neveu model [14]. In particular, the screening mass splittings can be solved from 
the reduced theory to an accuracy of g A {T), though only nonperturbatively. 



F. Reduced Lagrangian: bosonic part 

For the purpose of describing mesonic and baryonic screening states to order g (T), we 
only need the bosonic part of the Lagrangian Cb to order g 2 (T). In fact, these screening 
states reduce to two (three) valence quarks in infinite weak coupling (infinite T) limit and, 
therefore, the gluonic contribution to the binding coming from Lb has to involve one more 
loop and hence an additional g 2 {T) factor. For this reason, it is sufficient for our purpose 
to use the bosonic part of the reduced Lagrangian up to one loop, which has already been 



derived in Refs. 0,12 



Cb = ~\n v T% - \{D^)\D^f - l -m 2 D v> a + • • • , (43) 

where m 2 D = g 2 (T)(N/3 + Nf/6)T 2 is the one-loop Debye screening mass, //, v — 0, 1, 2 and 
the dots represent terms that contribute to the binding energy at orders higher than g A (T), 
e.g. the interaction terms between ip and A^. 



G. The complete reduced Lagrangian 

In summary, the complete reduced theory that should describe mesonic and baryonic 
screening physics with accuracy up to g A {T) for mass splittings and up to g 2 (T) for the 
overall mass is given by the sum of quark DR Lagrangian Lp in Eq. fl4"0|), the correspond- 
ing antiquark DR Lagrangian, which is still given by Eq. ([E]) with opposite sign for the 
couplings, and the gluonic Lagrangian £g in Eq. fl43|). The correspondence between the 
parameters in QCD and in the reduced theory is given explicitly in Eq. (|42|), and in Eqs. (§- 
|§p. In the kinematic region that is relevant for these screening states, i.e. those states that 
dominate the large- distance correlation between currents, the lightest quark mode is close 
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to its mass-shell, and the self-consistency of the reduced theory is guaranteed by the power 
counting rule given in Table |J, which have been justified in Refs. [16.17], exactly in parallel 
with the heavy quarkonium systems. The essence of the argument is that, if the interaction 
is weak, the off-shellness is small. 

A final word of caution about the usefulness of this reduced theory. Thermodynamic 
quantities, such as the free energy, specific heat and so on, are dominated by the gluonic 
zero modes. As it is well-known, quarks are heavy relative to gluonic zero modes and hence 
their contribution to thermodynamic quantities is suppressed at high T. This is perhaps 
the reason why quark degrees of freedom have never been considered relevant at high T, 
until one has been explicitly interested in the mesonic and baryonic screening states. The 
effective theory we derived is not intended for bulk thermodynamic observables, but only 
for those screening states that do not mix with purely gluonic states. 



H. At what temperature do we expect DR? 

At last we have obtained a reduced theory that includes the quark sector and that should 
be valid in the high-T limit. However, we still face the practical question of estimating the 
temperature above which this DR theory is going to be a good approximation to QCD. 

We know that DR is manifest only in a limited class of subtraction schemes [O], but 
even within this class there is some freedom left: one aspect of this freedom is the choice 
of the coupling constant in the reduced theory g\ = g 2 {T)T. While it is unambiguous 
how g 2 (T) runs with T, the numerical value of g 2 (T) at a specific temperature depends on 
how we match the reduced and the original theory. For instance, we have argued that a 
physically relevant way of choosing the coupling in the pure gluonic sector is to match the 
effective actions in the background field scheme ||19|| . This specific choice yields that DR in 
the pure glue sector sets in around 2T C , i.e. the coupling constant is sufficiently small at 
that temperature. 

This freedom in the choice of the relevant coupling constant can be rephrased in our 
formalism with the freedom of choosing the scale \i in Eq. (|7]) where c has different values 
in different schemes. Since our present interest is to solve screening states of quarks, it 
might be better to choose the subtraction parameter c with the criterion that the one-loop 
couplings #3 and g v be as close as possible to their tree-level values. For example, if we 
demand c\ — c 3 = 0, we get c = — 1/3 (if N = 3). This choice would make the subtraction 
scale n defined in Eq. larger than the one obtained from the criterion of optimal DR for 
gluons, which yields instead c = 1/22 (N = 3 and Nf = 0). Alternatively, if we demand 
ci — C4 = 0, we find a even larger subtraction scale // (c = —1.455). 

In either case, the fact that the subtraction scale is larger makes the effective coupling 
constant at a given temperature smaller in the quark sector than in the pure gluonic sector. 
This result is certainly in qualitative agreement with the empirical fact that the DR in the 
quark sector sets in at temperatures almost right above the critical point ||, whereas DR 
in the pure gluon sector sets in at about T « 2T C fl,|22"|. 
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VI. SUMMARY AND CONCLUSIONS 



We have shown with an explicit one-loop calculation that QCD at high temperature 
undergoes dimensional reduction also in the quark sector. More specifically, we have shown 
that static one-particle irreducible graphs contributing to the original four- dimensional cor- 
relation functions can be reproduced by a (2+l)-dimensional renormalizable Lagrangian to 
order g 2 (T) in the kinematic region where the lowest Matsubara quark modes are close to 
their "mass-shell" . Physical reasons why this kinematic region is relevant to screening phe- 
nomena have been also discussed. The reduced theory only contains the zero modes of the 
gauge field and the lowest quark modes as the explicit degrees of freedom and has the form 
given by Eq. (f|) plus the four-fermion term of Eq. fl20|) . 

Aiming at a better description of the mesonic and baryonic screening states, we have 
further improved the reduced theory to order g 4 (T) via a nonrelativistic reduction, which 
results in the reduced effective Lagrangian of Eqs. ([H]) and (^Bp. In fact, while the rela- 
tivistic version of the reduced theory mixes different orders in the coupling when used in 
the kinematic region relevant to screening physics, the nonrelativistic reduction explicitly 
separates the contribution of particles and antiparticles and allows a correct counting of the 
expansion parameter. 

Furthermore, we have also argued that the reduced theory in the quark sector, i.e. for 
screening mesonic and baryonic correlators that do not mix with entirely gluonic states, 
should become accurate at temperatures slightly above the chiral restoration transition 
temperature, due to the smallness of the appropriate running coupling Q 2 (T). In particu- 
lar, we find that the temperature above which the reduced theory becomes reliable in the 
quark sector should be even lower than the corresponding temperature in the pure gluonic 
sector [fT0[] . Our result has the potential for explaining present lattice data @,[| and provides 
a formal basis to the recent phenomenological modeling |6[] of the same data. 

We would like to stress that, although the reduced Lagrangian has been derived in a 
perturbative context, the reduced theory embodies all the infrared physics of the original 
theory, i.e. QCD at high temperature. Therefore, the solution of the reduced theory should 
reproduce the full long- wavelength screening physics of QCD in the high-T limit, which 
is nonperturbative: nonperturbative approaches such as lattice simulations are required to 
find this solution. Luckily since the large scale (T) has been explicitly factored out, it is 
now straightforward to put the nonrelativistic version of the reduced theory on a lattice 
following, for example, the method of Ref. |jl7|1 . 
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FIGURES 



FIG. 1. Vertices of C^ D that involve only light modes. A wiggly line represents a static gluon, 
while a thick (thin) solid line represents a quark of frequency uj + (w_), respectively. 

FIG. 2. Vertices of the original Lagrangian that are not present in £r D , since they involve at 
least one heavy mode (A£h)- A double line denotes a quark of frequency |c<j n | > 3ttT, while a 
spring-like line denotes a nonstatic gluon mode. The first vertex is "flavor-conserving" and the rest 
are all "flavor-changing". 

FIG. 3. Feynman graphs for the quark self-energy correction. Graph (a) is generated by ^C^ D , 
while graphs (b) and (c) involve "flavor-changing" vertices from A£h- The external quark carries 
four-momentum p = (cu + ,p) = (lo + , iuj + + q±, q 2 , (ft). 

FIG. 4. Feynman graphs for the quark-gluon vertex correction. Graphs (a) and (d) are gen- 
erated by vCpjj-), while graphs (b), (c), (e) and (f) involve "flavor-changing" vertices from A£h- 
The incoming quark carries four- momentum p = (uj + ,p) = (o>+, iuj + + qi, q2, Q3), arid the outgoing 
antiquark carries four- moment um p' = (uj + ,p') = (uj + , iuj + + q[, q' 2 , q' 3 ) , with \q\ <C T and \q'\ <C T. 

FIG. 5. Feynman graphs for the quark contributions to the vacuum polarization tensor. Graphs 
(a) and (b) are generated by £p>j3> while graph (c) involves "flavor-conserving" vertices from A£h- 
The external gluon carries four-momentum k = (0, k), with \k\ «T. 

FIG. 6. Feynman graphs for the one-loop correction to the composite operator. Graph (a) is 
generated by £r£), while graphs (b) and (c) involve "flavor-changing" vertices from A£h- The 
incoming quark carries the four-momentum p = (u;+,p) = (t<J+,iu;+ + qi,q2,Qs), and the outgo- 
ing antiquark carries four-momentum p' = (u + ,p') = (uj + ,—iuj + + q^q^q'^), with \q\ <g T and 
\q'\ <T. 

FIG. 7. Four-quark vertex (a) and its contribution to (b) one-loop correction of composite 
operators, (c) one-loop correction to the fundamental vertex and (d) one-loop correction to the 
quark self-energy. 
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TABLES 



TABLE I. One-loop coefficients. The entries that are not listed, i.e. those coefficients with at 
least one index equal to 2 or 3, are zero. 



coefficient analytic expression numerical value 

X - 24C'(-1) - (14/3) In 2 - 1 0.8898 

X 00 X + l 1.8898 

X 01 i[2 lE - 72C'(-1) - 10 In 2 - 6] i 0.1332 

X w i[X-l] -i 0.1102 

Xu iXoi -0.1332 

Y X 0.8898 

Too 2 7 £ + 24C , (-l) + (2/3)ln2 + 4 1.6465 

T i X 10 -i 0.1102 

Tio X w -i 0.1102 

Tn Xu -0.1332 

Z -2-/E + 12C'(-1) + (10/3) In 2 + 1/2 -0.3290 
Z 00 IE ~ 12C'(-1) - (7/3) In 2 0.9449 
Z i i[12C'(-l) + (4/3)ln2 + l] -i 0.0609 
Z w Zqi -i 0.0609 
Zn 7 g -36C / (-l) -51n2-3 0.0666 

Woo -2 7£ + 24C'(-l) + (17/3)ln2 + 21n3 + LT 1.1210 

W 01 -31n2 + 21n3 + LT 0.2383 

W u -24C'(-1) - (35/3)ln2 + 21n3 + 2 + LT 0.2011 



TABLE II. Power counting rules in the kinematic region appropriate for studying screening 
states at high T. 
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TABLE III. Coefficients of the one-loop corrections. 



coefficient analytic expression numerical value 

c[ C f [c - 3 + (X + In - iX w )/2] C f [c- 2.6768] 

c 2 C f [c-(3-X + X 00 + iX w )/2} C f [c - 2.0551] 

c 3 C f [c +(Y + Tn - iY w )/2] +C ad [c-1] C f [c + 0.3232] + C ad [c - 1] 

c 4 C/[c + {Y + Y m + iY m )/2] C f [c+ 1.3232] 



+C a(i [c - (T + Too + *T i)/4 + (Z + Z 00 + iZ i)/2] +C a(i [c - 0.3232] 
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